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Abstract 

The question of anomalous transport due to a band of impurity states in 
unconventional superconductors is discussed. In general, the bound state en- 
ergies are not in midgap, even in the unitarity limit. This implies that, gener- 
ically, the states associated with impurities are broad resonances, not true 
bound states. There is no impurity band in the usual sense of the phrase. 
The wavefunctions of these resonances possess interesting anisotropics in real 
space, but this does not result in anomalous hopping between impurities. 
I conclude that the system of resonances produces no qualitative modifica- 
tions to the T-matrix theory with impurity averaging which is normally used 
to treat the low-temperature transport of unconventional superconductors. 
However, users of this method often assume a density of states which is sym- 
metric around the chemical potential. This is not normally the case. It 
is found that the non-crossing approximation is not valid in a strictly two- 
dimensional system. 
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I. INTRODUCTION 



Impurity scattering plays a dominant role in the transport and thermodynamic properties 
of unconventional superconductors, far larger than in conventional s-wave superconductors. 
This is a consequence of the gap nodes which prevent the complete freezing out of scattering 
processes, and the fact that an anisotropic order parameter is far more sensitive to disorder. 
The critical temperature for an anisotropic superconductor is suppressed even in lowest 
order by the disorder potential. This follows from the breakdown of one of the conditions 
for Anderson's theorem ffl, which is that the momentum dependence of the pair potential 
is weak. A breakdown of the theorem leads to bound states in the gap when there are 
magnetic impurities in s-wave materials. It also leads to the possibility of such bound states 
from nonmagnetic impurities in the gap of an unconventional superconductor. This paper 
is devoted to questions about these states: their energy levels, their wavef unctions, their 
lifetimes, and the role (if any) they play in observable properties at low temperatures. 

This subject is topical because of the interest in high-temperature superconductors. Some 
of these systems appear to have gap nodes, implying the presence of an unconventional order 
parameter. No picture of these materials is complete without understanding the effects of 
dirt. Furthermore, experiments in the asymptotic low-temperature regime are special in 
that they probe the region of the Fermi surface near the nodes. It is in this region where the 
effects of impurities are most dramatic. The same considerations hold for the comparatively 
venerable heavy-fermion superconductors. Here we have solid grounds for supposing that 
some of these systems, particularly UPt3 and UBei3, are unconventional. Still, after more 
than a decade of investigation, the experimental details of the thermodynamic and transport 
properties of these systems at low temperatures are not fully reconciled with theory. In heavy 
fermion materials, however, it has become clear that strong impurity scattering, approaching 
the unitary limit, is the rule, not the exception. The Born approximation is inadequate. In 
high-T c systems, this is still under debate. In this paper, I will concentrate on this near- 
unitary limit. 
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The literature on impurity states in the superconducting gap begins with the papers of Yu 
and Shiba Q on magnetic impurities in s-wave systems. The two important ingredients 
are the pairbreaking nature of the disorder potential and the 'hard' energy gap - the density 
of states (DOS) is zero in some neighborhood of the chemical potential in the pure system. 
Bound states appear in the gap. Increasing the impurity concentration increases the number 
of bound states and decreases the gap, leading first to the gapless state and finally to the 
destruction of superconductivity {§. 

In unconventional superconductors, the bound states arising from ordinary potential 
scattering were first considered by Buchholtz and Zwicknagl ||. They concentrated on the 
Balian-Werthamer state, which has a hard energy gap, but the momentum-averaged gap 
vanishes: yjg A(fc) = 0, where the sum is over the Fermi surface. For such a gap, the results 
are somewhat similar to the previous case as the disorder potential is likewise pairbreaking. 
Although time-reversal symmetry is not broken, the randomization of momentum in the 
eigenstates prevents pairing by a momentum-dependent potential. Bound states appear in 
the gap. These authors also state that, in the unitarity limit of very strong potentials on 
the impurities, the bound states are at midgap. This statement has been repeated many 
times in the literature. However, I will argue below that it is incorrect. 

One important point about unconventional superconductors is that they cannot exist at 
very high impurity density. The critical temperature decreases as the impurity concentration 
is increased, and vanishes when %/r ~ ksT cQ where T c0 is the critical temperature in the 
absence of scattering. This implies that the regime of low impurity density is the only one 
of interest. 

With the discovery in the 1980's of the heavy fermion superconductors, there was an 
explosion of interest in the problem of disorder in unconventional superconductors. Many 
calculations of transport and thermodynamic properties at low temperatures have been 
published ||. The standard method, explained most completely by Hirschfeld et al. [[7], 
combines the T-matrix approximation with standard impurity averaging techniques. Gen- 
erally speaking, scattering near or at the unitarity limit is required to explain experiments 
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in both the heavy- fermion H, || and perhaps also in the high-T c materials 11, QTO | . This 



suggestion that the unitary limit is the appropriate one for Kondo lattice systems is due 



to Pethick and Pines ||11|| . While this model is certainly relevant to the weakly hybridizing 
/-level electrons in heavy fermion materials, its applicability to high-T c systems is unclear. 

The superconducting order parameters considered for both kinds of systems do satisfy 
t ne Eg = condition, but they do not have a hard energy gap. The DOS of the 

pure system is usually taken to vanish linearly or quadratically at the chemical potential 
/i. The standard method of treating the disorder potential leads to a finite DOS at \i |L2| . 



The neighborhood of the chemical potential where the density of states is flat is sometimes 
referred to as the 'impurity band' flOf . 

The impurity averaging method for unconventional superconductors has been explicitly 
questioned by some recent work |13| . The gap nodes lead to unusual wavefunctions for the 
bound states, with the possibility of anomalous overlaps between well-separated impurities. 
In compensated doped semiconductors, a high concentration of impurities can lead to a 
new conduction mechanism which predominates at low temperatures, conduction entirely 
through the impurity wavefunctions which form the impurity band []14[ . This possibility 
must be considered also in superconductors. The electrical conductivity of any such band 
would of course be shorted out by the conductivity of the condensate, but the opposite could 
well occur for the thermal conductivity. I will argue below that this does not occur. 

A more radical criticism of impurity averaging for two-dimensional sysyems is contained 



in papers of Neresyan et al. [[L5|], who find that multisite processes restore the vanishing of 
the DOS at the chemical potential. A recent preprint of Ziegler et al. shows that, for 
Lorentzian disorder, the finite DOS is not a consequence of impurity averaging. 

While the aim of the current work is to clarify the theoretical situation for nonmagnetic 
impurities in unconventional superconductors, there has been considerable recent work on 
magnetic impurities in both conventional and unconventional superconductors, stimulated 
by experiments [[IT]]. Some of this work has reached conclusions similar to those presented 



here, particularly with regard to the importance of carefully considering the real part of the 



Green's function in T-matrix calculations JT^|, [|19[ . 

In order to build up the theory from the start, I begin in Sec. [TI] with the question of 
bound states in the normal state of a semiconductor with a gap, in the limit of strong scat- 
tering. Since the semiconductor analogy is a powerful (but not omnipotent) one, this section 
provides much of the basis for the paper. The s-wave case is treated briefly in Sec. |m| , both 
to establish notation and to get a basis of comparison with unconventional superconduc- 
tivity. This latter topic, the main subject of the paper, is begun with calculations of the 
wavefunctions and lifetimes for single impurities in d-wave-type systems in Sec. |V[ Finally, 
in Sec. [V] the many-impurity case is discussed, along with the experimental implications for 
real systems. 



II. BOUND STATES IN SEMICONDUCTORS 
A. Introduction and formalism 

I examine an imaginary semiconductor in this section. The goal is to understand the 
process of binding an electron to an impurity with a very strong short-range potential. The 
physics of this process is sufficiently different from the textbook cases that certain features 
are likely to be as unfamiliar to the reader as they were to the writer. These features 
are important for the superconducting model which is believed to be of most relevance for 
high-T c and heavy fermion superconductors. 

Let us consider a semiconductor with a single impurity. The gap is the result of the 
lattice potential, a single-particle effect, and is not tied to the chemical potential, which 
lies in the gap. In the limit of weak scattering, this is essentially the familiar case of off- 
valence impurities in a Group IV material. This leads to impurity states very near the band 
edges. Our interest is in the opposite limit when the scattering is strong. The unperturbed 
Hamiltonian is 

ff. = Ew (i) 

k 
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The energies are measured from the chemical potential. I have omitted band and spin 
indices for clarity. The sum over momentum is always taken to include a sum over bands. 
The potential for a single short-range impurity is 



v = vJ2i^- (2) 

k,k' 

This is an s-wave potential. The phase shift is 

5(e) = -tan~ 1 (7riV (e)\/), (3) 

where N (e) is the density of states of the unperturbed system. The phase shift normally 
quoted in papers on transport properties is for states at the Fermi surface: 

5 = -tan- 1 [7riVo(e F )V], (4) 

The unitarity limit 5 — > tt/2 is reached when V — > — oo. 
The unperturbed Green's function is 

G Q (k,iu) = - — - — (5) 

The equation for the full Green's function for the Hamiltonian H = H + V is 

G(k, k',iu) = Go(k,iu)8^, + Go(k,iu)VGo(k', iu) + Go(k,iu)V^G (ki,iu)VG (k',iu) + ... 

fei 

= Go(k, iuj)5^p + Go(k, iuj)T(uj)Q (k', iu). (6) 

Here 

T(iu) = V + V 2 J2 Soil iu) + V 3 [£ Soft iu;)} 2 + ... (7) 

k k 

(8) 



1 - Vg (iu) 
and 



Continuing this function to the real axis: 

fc K fc K fc 

Thus 

Im g (u + i5) = -nJ25(u - eg) = -7rN (u), (11) 
where No(oo) is the density of states for one spin. Also 

Re g (u; + id) = V [ ^^W, (12) 
J cu — to' 

proportional to the Hilbert transform of the density of states. 
The expression 

T{u) = -, v ( \ (13) 
1 - Vg {u) 

shows that T has poles only when g is purely real and 

— = Re g {uj b ). (14) 

Let us agree that when real frequency arguments are used, a limit is implied where the real 
frequency axis is approached from above in the complex plane, corresponding to retarded 
functions. When Eq. [14] is satisfied but go has an imaginary part, then T is a Lorentzian near 
u b , and we are dealing with a resonance. If there is a pole, it represents a bound (V < 0) or 
an antibound (V > 0) state. In these cases we may write 

T-\u w u b ) w — - go(uj b ) -(u- uj b )g' (uj b ) = -(w - uj b )g' (uj b ), (15) 



where 



9oM = -Pj J — f p = -Z. (16) 

The integral for Z always converges because N (uj b ) = for a true bound state. Therefore 
we find 



in the neighborhood of the pole. The bound states are therefore characterized by poles in 
the T-matrix, and resonances by a sharp peak in the imaginary part of the T-matrix. 

The T-matrix gives the exact solution for the one impurity problem. It is not an approx- 
imation. 



B. Case of symmetric bands 

Let us consider a semiconductor with an unperturbed density of states which has a gap 
of width 2A, and satisfies the symmetry relation 

N (u)=N (-u). (18) 

I will argue below that this is not likely to be realized in the cases of interest, but it is the 
simplest mathematically. The density of states is illustrated in Fig. 0. Now we have 

N (uj')duj' 2oo 



Re g (co) = V \ 01 ' - « j, (19) 

J U! — U! UJg 



in the region |u;| << uj gi where 



N (u')du' 



< = J J > °' (20) 

uj g is of order A if the bandwidth is smaller than the gap energy (characteristic of insulators) 
and is of order the geometric mean of bandwidth times the gap energy in the other limit 
where the bandwidth is much greater than the gap energy (characteristic of semiconductors). 
The bound state energy ujb satisfies 

v = -^' (21) 

or 

2A 2 

u b = - — . (22) 
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For large \V\, (the unitary limit) this is a midgap state. This limit is shown in Fig. |l|. For 
V < the potential is attractive and the state sits just above the middle of the gap. This 
is an ordinary bound state. If V > 0, it sits just below the center of the band. It is an 
'antibound' state, but the wavefunction is localized, just as for a bound state. 



C. Density of states 

The density of states is 

N(u) = ImTr G(k,k', u + iS). (23) 

7T 

Comparing this with the equation 

g(k, k', a) = g (k, + g (k, ^)2»£ (£' , u ), (24) 

we find in the gap region (when Im Qq = 0): 

N(u) = N (uj) - i ~ ^k)' 2lm T (^) ; ( 25 ) 

k 

which may be written in terms of the change in the density of states: 

AN(u) = N(u) - N (u) = — - £kY 2lm T ( u )- ( 26 ) 

71 k 

Near uj = u>b this expression yields 

1 Z~~ l 
N imp = E(w - ^)- 2 }Im — = 6(u - u b ), (27) 

7T - U - UJb + 10 

k 

which is the impurity contribution to the density of states. To obtain the second equality, I 
have used Eq. |1^. In the region where Im go(u) ^ 0, we also have that Im T(u) ^ 0, and 
this represents a phase shift with an accompanying reduction of the density of states of the 
continuum such that 



AN(u))dw = 0. (28) 

The reduction of the density of states of the continuum just cancels the additional bound 
state (Levinson's theorem). 



D. Local density of states 



Near the bound state energy, the Green's function in real space has the form 



(29) 



u — u>b + iS 

Comparing with Eq. we may extract the wavefunction if)(r) by taking the Fourier trans- 
form of Gq{oj + iS): 

* • (30) 

uj b — e% + id 

Consider a semiconductor with a bound state at u>b and band edges at ±A/2. Let the bands 
be parabolic. Then the contribution from the upper band is: 

u h - (A/2 + kV2m) (31) 

poo pi gikrx 

= —Airm / k 2 dk / dx— -k (33) 

Jo J-i k 2 + fco 

4?rm r hdh elkr - e ~ ihr f 34 i 

~ ir Jo d k* + kl m 

= -^V^, (35) 
r 

where k^ = mA — 2mujb- The contribution from the lower band is the same except that ko is 
replaced by |mA + 2ma;&|. The wavefunction is very tightly bound, the decay length being 
short because the energy is far from the band edge. 

In the limit of a very weak attractive potential (V < and VN U << 1), then we are 
interested in the form of Re go(oo) when u « A. We find 

Reg (u) &N u log(— -)■ (36) 

The bound state energy is: 

u h = A - ce 1 ^. (37) 
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This is a state just below the upper band. The lower band has no effect in this case. It is 
important to note that the exponential dependence for the bound state energy is due to the 
fact that there is a finite jump in the density of states at the band edge. If there is a square 
root singularity: 

N (uj) ~ (cj - A) 1 / 2 (38) 

for uj > A, as one would expect in three dimensions, then there is a threshold coupling 
strength below which there is no bound state. 

The hydrogenic impurity case, of great practical importance, is different from all of these 
cases because of the long-range potential, which leads to an infinite number of bound states 
for all interaction strengths even in three dimensions. 



E. Case of asymmetric bands 

If the bands are asymmetric, N (u) ^ N (—uj), then the bound state is not in the middle 
of the gap even when \ V\ — > oo. This result is illustrated in Fig. ??. Consider an example 
where the lower band extends from — q to —A, and the upper band from A to e u . Let the 
bands have constant density of states and N u , respectively. Then 

ft S »w=»,r A ^ 7 «,r^ (39) 

J-e t UJ — UJ JAUJ — UJ' 

= -N e log | | + N u log | | . (40) 

If the enrgy is in the gap, < A, then this may be written as 

Re g (uj) = ^log|g| - iV tt log|^| - ^±^uj. (41) 
The bound state equation Re go(uj) = 1/V now has the solution: 



UJ b 



^ + AUog(|)-iV,log(|) 



(42) 



N u + N e 

Even in the limit of very strong scattering, this is not a midgap state. The bound state 
energy is displaced away from the band with the higher density of states because of level 
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repulsion. The asymptotic behavior of the wavefunction (the radius of the bound state) 
is still determined by the distance to the nearest band edge and the effective mass of that 
band. 

The effect of the asymmetry may be described as a renormalization of the potential in 
the following way. We may rewrite the eigenvalue equation as: 

uJb = (43) 

if we define the renormalized potential strength as 

V= V - s V . (44) 

H-VAf„log(i)-VAf ( log(J) \ + VN A y ' 

This equation defines the asymmetry factor Na, 

N A = N u log^£)-N e \og^£) (45) 

which is of the same order of magnitude as the density of states at the Fermi energy. It is 

V not V, that determines the energy of the bound state. It is important to note that, as 

V -> oo, 

V ±. (46) 

If the upper band is dominant, N A > 0, and the potential is repulsive, V > 0, then we have 
that V < V and the antibound state always stays below the center of the gap. This is simply 
a consequence of level repulsion. There is a similar effect for V < and < 0, with the 
bound state never reaching the center of the gap even if V — > oo. For the other combinations 
of signs, we will have a midgap state only in the 'accidental' case that V = —1/N A . 

This issue of band symmetry is crucial for the understanding of the bound state problem. 
It is particularly important to distinguish band symmetry from particle-hole symmetry, 
which is a very useful approximation for many calculations in superconductivity theory. 
Particle-hole symmetry is the assumption, approximately true in most cases, that the density 
of states of the normal material does not vary appreciably in the neighborhood of the Fermi 
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energy, the neighborhood being here defined as the range of energies within the cutoff energy 
huj c for the pairing interaction. The approximation may be stated as 



This is used in many elementary calculations of superconducting properties because only this 
range of energies is important for many purposes. A good example is the calculation of the 
critical temperature in the weak-coupling theory. The validity of the approximation arises 
ultimately from the mismatch of electronic and phononic (or other bosonic) time scales. 

Band symmetry is the assumption Nq(u) = Nq(-uj) which is essentially never valid. 
To give an idea of how far it fails, I have computed numerically the asymmetry factor for 
the following model semiconductor. It is a two-dimensional square lattice with a nearest- 
neighbor hopping matrix element t and a filling of 0.8 electrons per unit cell. The dispersion 
is eg = — 2.0t[cos(k x ) + cos(k y )}. At the Fermi energy ep = — 0.4t, there is a gap, symmetric 
around ep, of 0.021 The density of states is shown in Fig. |3|. Then Na is defined as 



where the brackets indicate an average over the band. The result is Na = 0.25/t per unit 
cell. Since the total band width is W = 8t, we see that the product NaW is of order unity. 
For any band, Na as a function of filling has one zero at some point. For this particular band, 
this occurs at half filling. In general, however, it is only for a special choice of ep that Na = 
and the band symmetry assumption is valid. The physical distinction between particle-hole 
and band symmetry is that there is no frequency mismatch for impurity scattering. The ionic 
potential which produces the band structure and the impurity potential are instantaneous. It 
is natural, but completely unjustifiable, to extend particle-hole symmetry to band symmetry. 
No conclusion which is based on such an extension is likely to apply to any real material. 




(47) 



N A =( 



1 



(48) 



13 



F. Level occupation 



In a semiconductor, the occupation of impurity levels is normally strongly dependent on 
the valence of the impurity relative to the valence of the constituent atoms. Here, we have 
been using a model in which the valence of the impurity and the background atoms is the 
same. We explicitly do not introduce additional states, only a potential which moves the 
old states around. This distinguishes the present work from Anderson magnetic impurity 
models, which generally do introduce such new states. 

The occupation at zero temperature is then as follows. If the potential is repulsive, the 
impurity 'peels off' one state from the valence band. It is therefore full, regardless of its 
position in the gap, and even if it is above midgap. If the potential is attractive, the impurity 
peels off one state from the valence band. It is therefore empty, also regardless of its position 
in the gap. In real semiconductors, it is normally true that the Coulomb repulsion prevents 
double occupancy of impurity levels, an effect not considered here. 



G. Many impurities 

1. Impurity band formation 

If there are N imp impurities at a finite density n imp in the system, we must consider the 
possibility that the wavefunctions on different impurities overlap. We begin with the case 
of two impurities. The Hamiltonian is then 

H = H Q + V l + V 2 . (49) 

Hq is the Hamiltonian of the pure system, V\ is the potential of the impurity at site r, and V2 
is the potential of the impurity at the site r + R. Our interest is in the limit Rk >> 1, where 
ko is the inverse of the bound state radius, as in Eq. [35[ The bound state wavefunctions 
satisfy 

(H + Vityi = u b ipi (50) 
14 



and 



(H + V 2 )if) 2 = u b ijj 2 - 



(51) 



The overlap matrix element is 



M 
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(52) 



The two terms are generally of the same size and asymptotic behavior. Taking the first as 
representative, we find 



M 12 ~ u b (il> 2 \il>i) ~uj b j vb*{r)vb{r + R)d\ 

where if) is the impurity wavefunction. Using Eq. |30], we find, for the symmetric case 

'e^d 3 k r e i% '< ?+ ^d 3 k'' 



(53) 



M 12 (R) ~uo b 



d r 



U) b 



Ak-R 



■ R d 3 k 



L0 b 



A'n 



(54) 



As in Eq. [35], we have the decay length kfi = mA — 2m,Ub- We may now write a Hamiltonian 
for the many-impurity case in the basis of the bound state wavefunctions at different sites. 
The resulting impurity bandwidth is of order u;&exp(— kon^ 3 ). u b is less than the gap 
energy A and l/fco is of the order of the lattice spacing. We generally expect a very small 
bandwidth for this 'deep impurity' (u b ~ A) case. There is therefore no metallic conduction 
when there is even a small amount of disorder in the impurity site energies. Interactions 
will also tend to localize the electrons and strengthen this conclusion. 



2. Impurity averaging 



The Hamiltonian for the many-impurity case is 



(55) 
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The impurities are located at the position The standard method of calculation is to 
average over the positions Ri (impurity averaging) |2(| . This method is valid for calculating 
the effects of impurities on the existing states if there are no correlations in the quantities 
ipQ(Ri)ipo{Rj) , where if>o(f) are the eigenfunctions of the Hamiltonian. This is a phase 
randomness assumption. The averaging process restores the translation invariance of the 
system on the average. The averaging method is clearly only appropriate when the number 
of impurities is an extensive quantity. 

It is convenient to rewrite the Hamiltonian in Eq. [55] as 

H = H + V + V - V , (56) 

where Vo is the spatial average of V. We then define Q Q (k,iu) as the unperturbed Green's 
function belonging to the Hamiltonian H + V . Both pieces of this Hamiltonian are diagonal 
in the momentum, and the second part gives only a rigid shift of the spectrum. The pertur- 
bation is then V — Vo, which scatters electrons from a state k to a state k'. The scattering 
amplitude is zero if k = k' because of the subtraction procedure. It is important to subtract 
the average potential explicitly, because the real part of the self-energy cannot be ignored 
in this problem, as it often can be in other contexts. 

The equation for the Green's function before averaging is 

Q(k, k', iu) = Q (k, ioo)8^p + Q (k, iu) V^pQ (k', iu) 

+ Q (k, iu) ^2 Vk,k»Go(k", iu))Vp, t pQo{k' , iu) + . . . 
k" 

(57) 

Here 

= V H ^' yAl ~ NimpVS^ . (58) 

i 

It is evident that in the extreme low density limit where the impurities have no influence on 
each other, this equation will reduce to the single impurity case. Quantities in Eq. [)7] are 
averaged using the prescription: 
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^= I A-Ujd 3 R l A, (59) 

i 

where L is the linear dimension of the system (henceforth taken to be unity) and N imp is 
the number of impurities. This leads to 



Q(k, k', iuS) = G imp (k } iu)Sk#- (60) 

A Green's function diagonal in the momentum Q imp (k,iu>) describes a state with uniform 
density. This shows that the averaging procedure washes out the density fluctuations which 
the impurities induce in the ground state (and other states) of the system. Bound states 
are an example of such density fluctuations. Thus, there is no possibility of bound states in 
this approximation. 

Carrying out the averaging and neglecting diagrams with crossed lines |2(J (an approxi- 
mation which will be discussed below) leads to the equation 

Q imp (k,iuj) = Q (k,iuj) + g (k,iuj)n imp V 2 ^2G imp (k" ' ,iu)Q imp {k,iu) 

k" 

+ g (k, iu)n imp V 3 Gim P (k", iuj) ^ Gi mp (k" ', iuj)Q imp (k, iu) + . . . (61) 

V k" 

The series is most conveniently summed by defining the the self-energy 

E(k,iuj) = g~ 1 {k,iu) - G^ p (k,iuj), (62) 
which leads to the equation 

E(k, iu) = n imp V 2 E^-^-Z(kM]-\ (63) 
1 — V T,j;>[ iu - e p — ^{k,iu)) 1 
It is then seen that E is a function of frequency alone for the short-range scattering potential. 

For that reason, the equation is algebraic, not integral. In the Born approximation, the 

denominator in this expression would be absent: 

S(fc, iu) = n imp V 2 YJ iuJ - e k> ~ iu} )V X ( 64 ) 
k' 

Eq. ^ contains two physical effects. The band edges are moved inwards because of level 
repulsion coming from the impurity potential. This directly affects the real part of the self- 
energy. The states are broadened because of the disorder which means that fc-states are no 
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longer eigenstates. This directly affects the imaginary part of the self-energy. Both of these 
effects tend to close the gap. 

To estimate the critical value of the potential strength at which this closure occurs, we 
must first specify the model a bit more precisely. Let us take the symmetric model of Fig. [I]. 
Then we need only determine when, at midgap ( uo = 0), the density of states first becomes 
finite. At real frequencies, let us separate the real and imaginary parts of the self-energy: 

E(w) = £'H +iS"M- (65) 

Since we have 

^)--^E (t ,_ £ ,_^ )2 + g2 , (66) 

we may simply increase the density ni mp and scattering strength V of the impurities, two 
quantities which occur only in the combination ni mp V 2 , and determine when S"(0) ^ 0. 
First consider the Born approximation. Symmetry dictates that S'(0) = 0. The imaginary 
part of Eq. |63] at zero frequency is then 

S // (0)=n rmp ^ 2 S"(0)/(0,S"(0)), (67) 

where the function / is defined as 

^"'"''^/ (.V+Vw - (68) 

and the integral runs only over the energies for which the unperturbed density of states is 
nonzero. Eq. |6^ always has the solution E"(0) = 0. It develops a second solution when 

' /(0,0). (69) 



f(0, E"(0)) is a positive, monotonically decreasing function of the nonnegative variable E"(0). 
This implies that the critical value of the disorder is 

(n imp V 2 ) c = A (70) 



This is when the inverse relaxation time corresponds to the gap energy, as might be expected 
on physical grounds. 

The consequences of band asymmetry are important for the density of states. Consider- 
ation of Eq. ^ shows that the real part of the self-energy is non-zero: 

S(^ = 0)^— log(^), (71) 

TIT 6l 

and that the derivative of the derivative of the imaginary part also does not vanish at midgap. 
Eq. B6] then implies that the density of states also has a nonzero slope at midgap. I have 



computed numerically the solution of Eqs. ^4] and |66] for the symmetric and asymmetric 
models. The results are shown in Fig. [|. It is seen that the minimum in the density of states 
shifts away from the middle of the gap. The chemical potential is given by a quite different 
equation than that for the gap minimum. It does not coincide with the minimum. 

The calculation for the full T-matrix equation, Eq. Q is only a little more complicated 
and will be omitted here. The result is that the threshold for complete closure of the gap 
is unchanged, but the spectral weight in the gap is larger for the same amount of disorder. 
The conclusions about the symmetry of the density of states in the gap are also unchanged. 



H. Physical Picture 

It should be evident that there are profound differences between the bound state cal- 
culations and the impurity-averaging calculations. The latter take into account only the 
potential fluctuations and even the sign of V is not very important. In the Born approx- 
imation, only V 2 enters the theory; even when the T-matrix is used, the main effect is to 
alter the extent of the phase randomization, not to create bound or antibound states. The 
scattering perturbs and broadens the extended states. The gap fills, but the states which 
are in the gap are extended states. There is no question of impurity band formation. 

In the bound state calculation, the effect of the impurity potential on the extended states 
is to give them a phase shift. This does not move the band edges. The bound states are 
split off from the bands. The gap fills with localized states. 
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Very strong potential scattering in systems with a band gap does lead to deep impurity 
levels. There is no reason for the energy of these levels to be at midgap. This could only 
occur accidentally. The decay of the impurity wavefunctions is very fast, as they are well 
split from the band states. This prevents the formation of impurity bands in all except 
pathological cases. 

How may we combine the results of both kinds of calculations, the bound state and the 
impurity averaging ? Let us first consider the set of all diagrams for the single-particle 

— » — » — » 

Green's function Q(k, k f , to) associated with the perturbation Hamiltonian J2i V(r — B4). If 
the impurities are numbered from 1 up to Nj mp , then a diagram of n-th order perturbation 
corresponds to an ordered sequence of these integers with n members. A given integer may 
appear more than once (repeated scattering from a single impurity). Intermediate momenta 
must be integrated over. Each scattering contributes a phase factor associated with the 
position of the impurity. The T-matrix calculations correspond to keeping sequences in 
which the same integer is repeated in succession many times. If the impurities do not 
interact, then only sequences containing a single integer are considered. This allows the 
bound state pole to form. The impurity averaging calculations omit the phase factor and 
omit the momentum integration after the final appearance of any integer. This is impurity 
averaging. The neglect of phase information prevents the build-up of the bound-state pole. A 
further approximation is to discard sequences in which numbers are interleaved, (an example 
is 1,2,1,2). This is the noncrossing approximation. It will be defined and discussed more 
carefully in connection with calculations in the superconductiong state. 

In the dilute limit, the T-matrix approach is appropriate if uj lies in the gap (and the 
unperturbed Green's function decays exponentially in real space), and impurity averaging 
is appropriate if u lies in one of the bands. It is therefore reasonable to conclude that the 
DOS is properly given by superposing the DOS from the two types of calculations, with the 
proviso that there are N — N imp broadened band states and N imp gap states. Here N is the 
total number of orbitals, proportional to the volume of the system. As the impurity density 
increases and the gap closes up, this energy separation argument no longer works. The 
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continuum states overlap in energy with the bound state. The phase information eventually 
becomes less important as the impurity states mix with the continuum states. The most 
likely scenario as the density is increased seems to be as follows. 

For low concentrations, the bound states will form a very narrow band in the gap. The 
integrated spectral weight in this band proportional to n imp . For subcritical disorder, there 
is no conduction associated with this band. There is some spectral weight in the gap because 
of the broadening of the band states. This weight is proportional to the total number of 
electrons, not the number of impurities. However, it does not overlap in energy with the 
bound state energy. If the disorder is above critical, then there is spectral weight everywhere 
in the gap, because of broadening and shifting of the band states. They overlap in energy 
with the bound state energy. This broadens the bound states, turning them into resonances. 
The qualitative behavior for n imp V 2 < {n imp V 2 ) c and n imp V 2 > {n imp V 2 ) c is shown in Figs. 
|5] and for the asymmetric case. 

The resonances can, in principle, play a role in conduction, since they are not necessarily 
localized. However, their density will normally be low compared to the band states, whose 
number is proportional to the total number of sites, not the number of impurities. This also 
means that if they do conduct, it is by hopping first into the continuum and perhaps later 
onto another impurity, not directly by impurity-impurity hopping as in an impurity band. 

The question of localization of all these states is subtle. In the impurity-averaging 
method, the single-particle Green's function is independent of position, seemingly indicat- 
ing extended states. However, when transport properties are calculated using this method, 
localization may appear in spite of this. Thus the states in the gap whose number of states 
is proportional to the total volume may or may not be Anderson-localized by the disorder 
- the impurity-averaging calculation of the single-particle Green's function as carried out 
here gives no useful information about this. Conventional wisdom tells us that, in the low 
impurity concentration regime considered here, states near the chemical potential should be 
extended in three dimensions and localized, but with extremely long localization lengths, in 
two dimensions. 
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III. S-WAVE SUPERCONDUCTORS 



The calculations of the effects of impurities in superconductors are very analogous to the 
calculations in semiconductors. This analogy is most easily exploited if we introduce the 
Nambu operators. The defining equation is 

The Pauli matrices t±, T2, t 3 and the identity matrix t (omitted when clarity requires), act in 
this two-dimensional space. The unperturbed mean-field Hamiltonian for a superconductor 
with a constant gap A is: 



= E *hr-s% + E *I-ATitt s , (73) 
k k 

and the impurity potential is: 

k,k> 

The matrix Green's function is defined as 

g^iu) = -<T T ^T)^(0)>. (75) 
In comparison to ordinary notation we find 

— # — # 

C?ii (fc, iu;) = £ TT (A;, iu;) (76) 

022 (k, iw) = -Gu 0, -iu) (77) 

Q 12 (k,iuj) = T{k,iuo). (78) 

Since the density of states is independent of the spin direction in singlet superconductors, 
we have that 

N(u) = --Im Tr Gu(k,u + i5). (79) 

7T 

We may now easily calculate the unperturbed Green's function: 
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Go(k,iu) = - — 

iu - e^r 3 - Ati 

The T-matrix is again denned by: 



iu + ep- 3 + Ari 
cu 2 + e 2 r + A 2 ' 



(80) 



G(k, k', iu) = G (k, iu)5^ p + Q (k, iu)T(iu)G (k' , iu), 



(81) 



which leads to: 



T(iu) = Vt 3 + V 2 r 3 g (iu)r 3 + V 3 T 3 g (iu)r 3 g (iu)T 3 + 



(82) 



where 



c \ ST r (Z ■ \ V" iuJT ° + e k T3 + Ari 

0(>M = L IW) = - 1, ^2 + e 2 + A 2 • 

fc k k 



We shall assume particle-hole symmetry so that 



H = 0, 



(83) 



(84) 



where u c is the cutoff frequency of the interaction, which is much smaller than the bandwidth. 
We also recall that in weak coupling A = for \e%\ > u c . Then: 

de 



ru 

go{iu) = -{iuT Q + Ar 1 )N (e F ) / 

J — t 



u 2 + e 2 + A 2 



(iuT + er 3 )N^e)de 
u 2 + e 2 



i\e\>u) c 

-{iur + Ari)iVo(eF 



de 



u c u 2 + e 2 + A 2 



+ A s 



= -(lUT + An ; = V == + A, 



Vco 2 + A 2 



(85) 
(86) 
(87) 



The asymmetry factor is defined as 



A s = T 3 N A = T 3 



E 

L k 



- 2 + 4 



(88) 



where the sum runs over wavevectors whose energies are farther from the Fermi surface than 
the cutoff frequency u c . We find 



AT A = AUog(^)-iV,log(% 

u c u c 



(89) 
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in the model of constant density of states beyond cu c . A s also has an imaginary part, but it 
is smaller by a factor of u/u c . 

We may now compute the T-matrix: 
T (iw) = Vr 3 + V 2 t 3 (clt + bn + N a t 3 )t 3 + V 3 t 3 (clt + b n + N A r 3 )r 3 (ar + br x + N A r 3 )r 3 + . . . , 

(90) 

where the definitions a(iu) = — ircNQ(eF)uo / \/uj 2 + A 2 , and b(iu>) = — ttN (€f) A/\/uj 2 + A 2 
have been made. It remains to sum the geometric series and invert the resulting matrix: 

oo 

T(iu) = Vr 3 J2[V(ar + br x + N A r 3 )r 3 ] n 

n=0 

= Vt 3 [t - Var 3 - Vbr x r 3 - VN a t^ x 

= Vt 3 [t - Var 3 - Vbnr^ 1 

aV 2 T - bV 2 n + t 3 V 
~ 1 - a 2 V 2 + b 2 V 2 

In these equations, we have made the definition, as before: 



(91) 



V = . (92) 

1 + N A V 1 J 

This shows that the binding potential is renormalized by the band asymmetry just as in the 
semiconductor case. This limits its strength to something on the order of the bandwidth. 
Continuing to the real axis, we find 

= a(w)V 2 T - 6(u)V ,2 7j + t 3 V 
y ' 1 - a{uo) 2 V 2 + b{oo) 2 V 2 V ; 

ajuJ^ro-biu^n + r.V 
l + (7rN (e F )V) 2 

This function is nonsingular, indicating that there are no bound states in the gap. This 
arises from a cancellation of the frequency dependence in the denominator, essentially that 
of Anderson's theorem. 

The asymptotic dependence of the Green's functions at large distances is of interest for 
later considerations. The computational method is very similar for all cases, so it is given 
in some detail here and abbreviated later. In three dimensions, we have: 
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G(f,iu) = I e lk - r g{k,iu)d 3 k (95) 
No f de f 1 d X d ^ + ^ + ^ e*". (96) 



4vr 7 7-i Jo cj 2 + e 2 , + A 2 . 

k 

The x- integral is performed by choosing a contour which runs from z = — 1 + 200 to — 1, 
then from —1 to +1, then from +1 to 1 + zoo. Since we will use this contour several times, 
it is shown in Fig. [F] for reference. Let u = (&p + e/vp)r. This variable must be treated 
carefully as the approximation of linear dispersion is a very limited one. We may extend the 
limits of integration over u only if the integral is rapidly convergent. This usually means 
that u must be integrated last. Bearing this in mind we have: 

/ dze iuz = f 1 dxe iux + i r dye iu{1+iy) - % H dye iu{ ~ l+iy \ (97) 
J J-i Jo Jo 

as long as u > 0. The contour integral vanishes because the function is analytic, so we find 



/l poo /*oo 

dxe iux = -ie iu / e~ uy dy + i / e~ uy dy 
-l Jo Jo 



--(e iu -e- lu )^--^(e iu -e- tu ). (98) 

U KfT 



Substituting yields 



Q{rM = [ de[e^M _ e -ir(k F+ e M] ^ + J + An 
2k F r J uj l + e z + A 2 

mN e- dr ' VF 



2kprd 



-[(iujTo + idr 3 + An)e lkFr - (iuT - idr 3 + Ar 1 )e" l ' CFr ]- (") 



Here d = \/ lo 2 + A 2 . As A — > 0, this becomes 

Q(f, iiv, A = 0) = -lj^e-M r ' v *[{sgn{u)T + r 3 )e^ - (sgn(co)r - T 3 )e- ik ' r ]. (100) 

This is the normal state Green's function. The off-diagonal components are zero, and the 
diagonal one may be written as: 

g u (f, iu, A = 0) = - 7 L^l e -\^hF e ik F r SS nM_ ( 10 i) 

kpr 

For future reference, the corresponding result in two dimensions is 

g il (f,itU,A = 0) = e -Wr/v Fe ik F r sgn(u>)_ ( 1Q2 ) 

\j2ixkpr 
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At low frequencies, u « A, we find 
g(f,iu,u« A) = -^L e - Ar ' v >>[(sgn (w)n + r 3 )e^ r - (sgn (w)n - r 3 ) e -^]- (103) 
The particle-hole part in the low frequency limit is 



g 11 (?,u,u« A) = -vrA^/-^^. 



(104) 



The density of states is zero in the gap, so the Greeen's function is purely real. The 
exponential damping in real space is also due to this fact, the pole of the function 

1 



e 2 + u 2 + A 2 

being off the real axis. The decay length is the coherence length i>f/A. 



(105) 



IV. UNCONVENTIONAL SUPERCONDUCTORS 
A. Bound state energies 
I will restrict the discussion to the singlet case. The new feature in the equations is that 

— * 

the gap function is now /c-dependent, and satisfies: 



k 

The Green's function is the same as above, and the T-matrix is still given by 



(106) 



T(iu) = Vt 3 + V 2 T 3 g (iu)r 3 + V 3 r 3 g (iuj)T 3 g (iu)T 3 + 



(107) 



where 



k 

= -iuN r ^- I dVL k 



= -iujT 7rN — 

4:71 



de 



OJc 



cu 2 + e 2 + A| 



de 



\e\>u> c uj 2 + e 2 



V^TAf 



(108) 
(109) 
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Here fig is the solid angle on the sphere. The band asymmetry is represented by the second 
term in the equation, and comes from the integral over energies far from the Fermi surface. 
Its effect is the same as above, i. e. , it results in the replacement of the bare potential V by 
V = V/(1 + N A V). 

The equations for the unconventional case are, to this point, actually somewhat simpler 



than the s-wave case. The T-matrix is obtained from Eq. ^Tjby setting 6 = and a = go(iuj) : 



1 - g (iu) 2 V 2 



Component by component, we find 



TuM = ^Yl±l = yl±Sm =V 1 , (111) 
l- 9 oM 2 f ! 1 -V'nU^) 

for the up-spin electrons, and 

r B N = v ( 112 ) 

for the down spins. The corresponding equations for the binding energies are: 

Reg {u b ) = ^, (113) 

and 

Re g (-u b ) = -I. (114) 



Since, from Eq. |108| , Re go(uj) is odd (if we exclude the asymmetry term), these equations 



heave the same solution for the bound state energy uj b , as indeed they must. 

To get some insight into the equations, we compute the T-matrix for the 'polar' case, 
which has a line of nodes on the equatorial plane of the Fermi surface: = A cos(6 1 ^). We 
have 

iVn f°° r 

Im g (co) = -tt-^ / de / dQ^[S(co - E % ) + 5(cu + E % )] (115) 

47T J— oo J 

= -vriVo [°° de [ dx6(U - Je 2 + A 2 x 2 ), (116) 

J-oo JO 
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which gives 



ImU<»)= 'o'f 07 ' 2 ' M <A o ( 117 ) 

ZZA 

sin ( — ), > A . (118) 



A LJ 



Computing the real part: 



Re 9 „M = -\J Im9 J { ^' ("9) 
2A V-a — a;' 2 ja — a;' J-u c lu — lu' 

log I ^ o ^ log 1 7 ; r-r f| (121) 



2A "'Ag-cu 2 ' 2 & 1 (w c + lu)(A - lu) 

N TTU ( 1 

AT ( 2 g 



a; 2 



(122) 



at small w. The asymmetry term A s has not been explicitly included in this expression, 
since we may more conveniently include it in the potential strength. [A s = —N Q \og(e u /ee) 
in the simplest model density of states.] For other gap functions which give an even smaller 
density of states at low energies, the slope of the real part is completely dominated by the 
first term, and gives a generic result: 

Re g (u) ~ -P^, (123) 

where (3 is a number of order one. The bound state energy is: 

oJ b = (124) 



as long as NqV » 1. This is the more generic result from Eq. |123 , which we will also use 



below. The logarithmic corrections in the polar case are only important if, for accidental 
reasons, 1 + N^V is exponentially small. 

Once again, in this case, the bound state energy does not go to zero even in the unitarity 
limit. This is very important for the properties of the bound state wavef unctions, to which 
we turn next. 
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B. Asymptotic Spatial Dependence 



At large distances from the impurity, the wavefunctions for bound states are determined 
by the asymptotics of the diagonal elements of the real-space Green's functions evauated at 
Ub [30[ The off-diagonal elements are also of interest, as they represent the suppression of 
the gap in the neighborhood of the impurity. Their asymptotic spatial dependence is always 
the same as that of the diagonal elements, so I do not give expressions for them separately. 

In three dimensions I shall concentrate on the polar case. The gap takes the form 
A(k) = A cos#£. We shall take a spherical Fermi surface; for this situation the gap has 
a line of nodes around the equator. As in the previous section, we assume that the pair 
wavefunction is a singlet in spin space. This total pair wavefunction does not satisfy the 
Pauli principle. However, it is the simplest example of an unconventional state with a line 
of nodes in three dimensions and the calculations for it are already complicated. 

Now take the direction f from the impurity to be in the z-direction. I will work out the 
Green's function in the main text for this case. For other cases, the intermediate steps are 
relegated to the appendix. 

Then 

The important integral (because it expresses the angular dependence) is the x integral. 
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This is: 



I(e,iu b ) = J 



W 2 + e 2 + A 2 x 2 flX 
A i-i fi 2 + s 2 + x 2 

A J e Q2 + s 2 + z 2 aZ 

_. iu J_ r-u„ M + STs + (l + iy)n 
A Jo n 2 + s 2 + (l + n/) 2 y 

+l€ A Jo Q 2 + s 2 + (-1 + iy) 2 V - 1 ' 

The method is the same as in the s-wave case, and we use the same contour, as shown in Fig. 
|7]. The scaled variables s = e/A , Q = u^/Ao, have been used, and again u = {k F + e/v F )r. 
We are interested in the limits Q << 1, and k F r » 1. The contour integral is easily 
performed by the residue theorem, and the strongly peaked function e~ uy allows us to do 
the other two integrals. Thus, 



^ + sr 3 + iy/Q 2 + s 2 Ti . e m sr 3 + t\ . e IU sr 3 — t\ 



I(e, - Aq e u 'uA 1 + s 2 + fi 2 + \a 1 + s 2 + fi 2 ' 

(127) 

For u >> 1, we have the useful integrals: 

^Lw^s 1 '^' (128) 



2tt|0| u 

-oo ' li 



h = r e~ u ^ + ^ 1/2 ds » (ffl)iA» e -«mi (129) 



This yields 



/•oo 

/ ds/(e, u; 6 ) = — e 

J — oo Z\n 



u|n| 



^(4^) 1/2 + ^(^) 1/2 

U Si It 



tt(S1 + VI + Sl 2 r 3 - in) 



+ V T y 7; ' 1J exp[tk F r - VTTIV 2 A r/v F } 
A^kpry 1 + ii 

+ ^- a ;f n ^ElLT' T jl eM -ik F r - VT+^A r/v F ]. (130) 
A k F rV 1 + S2^ 



The Green's function at large distances is therefore 
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Q(rz,iuj h ) = — / dsl(e,u b 

2 J — oo 



2 A kpr 

7rJV e-K a +^) 1/a ^ 
2(u; 6 2 + A§)V2 X 

[(w 6 + K 2 + Ao) 1/2 r 3 - in)e ik - r + (-co b + (co 2 + A 2 ) 1 / 2 ^ _ lTl )e'^]. (131) 
The component of most interest is the particle-hole part at real frequencies: 



2 AqKft 

2 

[ ^ + (A 2 -c 2 )i/ 2 _ ^-(Ag-o;?) 1 / 2 \ 

\ (A 2 -c 6 2 )V 2 (A 2 -c 2 )V 2 /■ ^ 

The first term is the pole contribution, which arises from interference of the low energy 
states near the line of nodes. If the state is precisely at the chemical potential, then these 
states interfere to give a power law decay. There is no oscillatory component because the 

— * 

states on the nodal line have k ■ f — for this direction of r. The length scale of this decay 
is A /\ub\kF, which is the interatomic spacing unless the binding energy happens to be very 
small. 

The second term is the stationary-phase contribution which is due to the states near the 
north and south poles of the Fermi surface. Since these are gapped, we get an exponential 
falloff which reflects the energy gap at these points. There is an oscillatory behavior because 

— * 

k ■ r = ±/cpr at these points. The length scale of the decay is the 'local' coherence length 
(coherence length at the pole) A{kz)/vF if the binding energy is small. 

A second special case of interest is when the direction of r is on the plane of nodes. Then 
we have 

The details of the integration are given in the appendix. The result is: 

G(rx,iu b ) = Z^ e -I^|r/*F [jXr^sgn ( Ub ) TQ + r 3 ) - e - ikFr (sgn (cu b )r - r 3 )] . (134) 
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The particle-hole part is 

g U (rX,U b ) = Z^ e -H\r/v Fe ik F r sgn K)_ Qgg) 
2kpT 

These are essentially the normal state results. The point here is that the stationary phase 
points at k = ±kpx have no gap and the behavior is therefore entirely normal. 

The general case is rather complicated. Let the direction from the impurity be inclined 
at an angle ip to the z-axis. Let us define 

G(r,ip,iub) = Q(r (sin ipx + cos ipz,iuJb)) (136) 

We have that 

g{r - c -'^ T,L ^ + 6H %x-M-xT ' (137) 

where a = cosip, b = cos^sin and u = (kp + e/vp)r. The coordinates for the integration 
have been rotated so that the polar axis defined by x = ±1 is along r. The particle-hole 
part of the Green's function is found in the appendix to be: 

kpr \AoCos z ^sm ip I 



iN 0K I ik F r Ub+ V A c OS 2 ^-^ b 2 _ _ ikpr U b - V A c OS 2 V ~ 

,2 



2k Fr \ ^JA 2 cos 2 i) - uj\ ^A 2 cos 2 ^-w fe 2 

xe V A ? cos2 H r / , 'J', (138) 



This general form does not reduce to Eq. |135| when xjj = tt/2 or to Eq. |132| when ip = 0. It 



is only valid in the intermediate regime of angles. The result does, however, show that the 
two distinct decay behaviors noted above are both present in the generic case. 

C. Two dimensions 

We now turn to the case most relevant to high-T c , which is two dimensions with a singlet 
d-wave gap, (which is consistent with the Pauli principle). The Fermi surface is a circle and 
the gap will be taken as 
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A(k) = A cos2</>g. 



(139) 



This d-wave gap has nodes at the intersection of the lines k x = ±k y with the Fermi surface. 
Let us first determine the Green's function for the direction along one of the axes. We have 



g(rx,iu b ) 



2N 



de 



1 dv iu b + er 3 + A Ti(2v 2 — 1) 



(140) 



7T J J_ X v/l"^ w 2 + € 2 + A 2( 2v 2 _ 1)2 

Here we have defined v = cos</>g so that A = A (2t> 2 — 1). Also u = (hp + e/vp)r. The 
integrations are performed in the appendix. The particle-hole part is 
2iN 



Gu(rx,w b ) 



x 



7T 



7T 



i sgn [u b ) cos 



(Ag-u; 6 2 )V4 
k F r 



w " M/2 I 

)(— ) 1 exp 



kpr' 



Al-uol-AA x 



An - loI 



A f 



'1 - i)n 



+ 



2 V k F r 
'1 - i)ir 



W e ik F r e ~^/A2-u 2 b r/v F 



iQ + i\ Al 



A 2 , - ujI 
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-ik F r -y/ Al-uj^r/vp 



iVt — i\ A, 



10 



2 ~\ 



A 2 / > 2 



■ 
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2 y A; F r 

The peculiar phase factors reflect the orientation of the nodal directions relative to the 



crystal axes. The overall behavior is very similar to the three-dimensional case, Eq. |1 32 
The direction along the node is more easily computed. 

2N 



Q\J — —,iu b ) 



de 



dv iu b + er 3 + A TiV\/T 



v ■ 

,2\ ' 



(142) 



V2 ' °' vr J J-i cu 2 + e 2 + 2A^ 2 (l-t; 2 

where the part proportional to T\ vanishes by symmetry. The integrations are performed in 
the appendix, with the result that: 
2iN 



r i x+ y \ 

Qu{r—^-,uj b ) 



x 



7T 



4W2 , sgn M ( / A TT3 _ A y/^ e -k F r { ^^ b -^ /V ^ 



k F r ' (A 2 -^ 2 )V4W o 
^-SB. ( sgn (a*) + l) e ~ 

+ (1+Ott £^ ( ^ ( ^ _ 1)e _|. 6 | r/ „ F 



y/k F r 



(143) 
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This is essentially the two-dimensional normal state result. 

Now consider the general case, where the direction is inclined at an angle ip from the 
diagonal and define: 



7T 7T 

Q(r, ip, iu b ) = Q(iu b , r cos(- - tp)x + sin(- - ip)y). 



(144) 



The result, derived in the appendix, is: 



r ( i \ 2iN ° 
Gii{r,ip,u b ) = x 



7T 



8tt 3 ^2 



1/2 sgn u b 



y k F r\ cos V| (^o-^ 2 ) 1/4 

-fc F r|co S ,/.|( v /A2- w 2-A () ) 1 /2/ v ^A^ 



I7I (V / A^-Ao) 1 /^ 4 x 



cos 



tk F r sinMAl - - Ao^/V 2 ^ 



(1 — j\qr p ik F rcosip 



2 \Jk~Fr\ cos VI 



f 1 ±.'W „-ik F rcosip 

(1+zjTT e ( ^ ^ _ 1)eH . bCOS ^| r/ „ F] y 

fcpr| cosVl 
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D. Discussion 

As compared with the semiconductor, there are several interesting differences in the 
bound state wavefunctions of the unconventional superconductor. The most important is 
the fact that the decay length may become longer, not shorter, as the bound state energy 
approaches midgap. In both cases, however, there is nothing special about the middle of 
the gap in a real system, so exponential decay is still the norm. 

In the superconductor, the wavefunctions have two components when viewed in real 
space. The two components correspond to two different exponential decay lengths. One of 
these lengths is determined by the gap along the direction of propagation. This length is 
Vpj \J A(k) 2 — ujI. This is similar to the s-wave case except for the anisotropy. Indeed, this 
contribution comes from the fully gapped region of the Fermi surface. The second length is 
VF/\u b cos VI; where is the angle away from the nodal plane. This contribution is peculiar 
to the unconventional case, arising from the gap nodes. These lengths are anisotropic, so 
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the wavefunctions are also anisotropic, with 'arms' in the directions of the gap nodes. Some 



very nice pictures of these wavefunctions may be found in Ref. [jTS . 

The decay is always exponential in all directions unless, for accidental reasons, the bound 
state energy is zero. 

V. MANY IMPURITIES IN SUPERCONDUCTORS 
A. Introduction 

The discussion of many impurities in the metal and the superconductor are usually 
considered to be parallel, and the same equations, with only the generalization to the Nambu 



formalism, are used for both [20|. In articular, the method of impurity averaging is not 
modified. Hence, in this section we shall not repeat the calculations of Sec. O G| . 

For unconventional superconductors, without a hard gap, this is fundamentally reason- 
able. We have seen that impurity averaging becomes valid when the continuum states 
overlap the bound state energy. The constructive buildup of phase required to make the 
bound state is destroyed when there is overlap with the continuum states, themselves pos- 
sessing a random shift. This is always the CcLSC, clS the bound state energy is not in midgap. 
Accordingly, the pole found in the T-matrix calculations should immediately broaden into 
a resonance. This being said, one should address four basic issues which arise in practical 
calculations. These are: symmetry of the DOS around the chemical potential, the nature 
of the states at the chemical potential, the validity of the noncrossing approximation, and 
impurity band formation. 

The first two are relatively easily dealt with. The third and fourth are treated in the 
next two sections. 

The density of states is not symmetric about the chemical potential, even when particle- 
hole symmetry is valid. This has already been pointed recently by other authors in the 



magnetic impurity case [18], [19]. This arises from the same source as in the semiconductor. 
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Symmetry of the DOS requires that the bands themselves be symmetric over their whole 
energy range, not just over the neighborhood of the chemical potential. It appears that this 



fact is often not taken into account in practical calculations. The classic reference |20| advises 
us to neglect the real part of the energy shift. This is indeed safe for the neighborhood of 
the Fermi energy in a metal, and for weak scattering in s-wave superconductors, the cases 
discussed in Ref. pOfl . It is not valid in unconventional superconductors with a soft gap. 

The states at the chemical potential are sometimes termed 'bound states', and their 
heritage as the descendants of the T-matrix poles is emphasized. It should be clear from the 
discussion that this is not correct. These states are the broadened and shifted continuum 
states. The number of such states in any range of energies is proportional to the total 
number of orbitals in the system. The daughters of the bound states will generally live in 
a resonance away from the chemical potential. Their number is proportional to the number 
of impurities. 



B. Noncrossing approximation 

One approximation used in nearly all calculations of superconducting properties is the 
noncrossing approximation. This has been questioned in recent work, ||15|| , and I reproduce 
and expand this criticism here. This approximation is defined diagramatically by repre- 
senting each impurity as a cross through which momentum flows and is conserved. Let us 
restrict the discussion in this section to the Born approximation for simplicity. Then two 
typical diagrams for the normal state are shown in Fig. || We take a circular Fermi curve in 
a two-dimensional system. Diagram (a) in Fig. || has no crossed lines, whereas diagram (b) 
does. The same processes are shown in momentum space in the diagrams in Fig. |9|. Diagram 
(a) describes a retraceable path, while diagram (b) contains a circuit. Therefore, the sec- 
ond diagram must satisfy one additional momentum conservation condition. It is therefore 
smaller in magnitude than the first. Explicit calculation shows that the small parameter 
involved is (l/kp£), where t is the mean free path. 
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In a two-dimensional d-wave superconducting state at low temperatures, this argument 
must be reconsidered. All diagrams must effectively satisfy additional momentum con- 
straints, as only scattering between gap nodes, situated at (±fcp/\/2, ±&f/v^2) is important. 



The two diagrams in Fig. [10] again correspond to the processes of Fig. |8]. The two diagrams 
are of roughly equal weight, even though the second one is crossed. There are no additional 
constraints which must be satisfied by the second diagram. 

The noncrossing approximation is therefore very questionable in unconventional super- 
conductors with point nodes. The authors of Ref. ]15| attempt to go beyond this approx- 
imation in a not very realistic model. No calculations of transport properties have been 
carried out except using the noncrossing approximation. It certainly is difficult to justify for 
the d-wave states considered in the context of high-temperature superconductivity, if these 
materials are taken to be two-dimensional. 

C. Impurity band formation 

Does there exist the possibility of the formation of an impurity band ? Does conduc- 
tion in this band influence, or even dominate, the transport properties in the limit of low 
temperatures ? 

We begin, as in Sec. [II G| , by considering two impurities. The overlap now depends, to 



some extent, on the direction of the vector connecting the two impurities, with the direction 
of minimum gap being the direction of maximum overlap. If we take the two-dimensional 
example summarized in Eq. |145| the overlap proceeds according to exp(— kpr\ cos^l) or 



exp(— \u b cos if) | r/vp). 

The introduction of many impurities always brings one new number to the problem: the 
average distance between impurities, which shall be denoted by £ imp ~ where n imp 

is the number of impurities per unit volume. If the impurities are not identical, we have 
a disorder parameter W, defined as the width of the distribution of the potential strength, 
previously the single number V, of the impurities. The usual model of impurities is that they 
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are all identical: W = 0, but we will consider also W 7^ 0. If interactions on the impurities 
are important, we may introduce a Hubbard-type parameter to describe the interaction 
strength. We shall not discuss this possibility, but only note that this also introduces a 
breaking of the band symmetry which can move the bound state away from midgap. This 
situation is treated in Ref. ITCH . 



The first question for band formation is the following: given a wavepacket located at an 
impurity site, is it more likely to hop to a neighboring impurity, or to leak into a continuum 
state ? If the latter, then the impurities merely form a system of resonances and averaging 
procedures should be approximately valid. In this case, we may make arguments similar to 
those for the semiconductor to argue that the results of the two types of calculations may 
be combined. 

If W = 0, then the bound state energy ujb is fixed at some position relative to the chemical 
potential fi. In the general case Ub 7^ fi, even in the unitary limit, as we have seen above. 
This means that there is a finite density of states at u>b- The lifetime of the wavepacket for 
decay into the continuum r c is given by 

— = 7riV(^)|Vf. (146) 



T, 



c 



Near the unitary limit, we have \V\ > 1/N . Using Eq. |123| , we then find 



1 &h / \ 

This may be anomalous only if \ujb\ « A . The rate for interimpurity hopping 1/Tj, is of 
order 



(148) 



n N 

where £<j is the minimum decay length. As shown above, we have that ^ = vp/\oJb\- 
The criterion for band formation is then 

e -Wfd> Z^. (149) 
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Unless u>b is accidentally very small, this means that we must have £i mp ~ However, this 
is the dirty limit. This limit does not exist for unconventional superconductors because the 
critical temperature T c is a sensitive function of impurity concentration and the situation 
d-imp ~ id corresponds to T c — ► 0. Hence band formation does not occur in the W = case. 
The only possible exception would be if, by some accident, \ujb\ << A. This is an intriguing 
possibility. However, it is not related to the unitary limit. 

Now consider the case of finite W. There will be a distribution of bound state energies. If 
this distribution does not include the chemical potential, then the previous conclusion that 
no impurity band forms remains valid, as no quasi-bound states have low enough energy to 
be anything but broad resonances. The interesting case is when the distribution is broad 
enough that some of the bound states have very small \ujb\, (|^&| << A). 

We may build up the state by considering pairs of impurities. Almost all such pairs 
which involve a low-energy impurity state (energy |u;&i|) will then involve as the second 
member a state for which \u\a\ ~ A. These states will mix, with overlap matrix element 
M 12 . Under the influence of the mixing, state 1 then has probability amplitude on site 2 of 
Mf 2 /(ubi — u;^) 2 , when |M 12 | << \{oJbi ~ UJ b2) I- The transition rate to the continuum for 



this state is then 

1 M? 2 tt/3 



(150) 



t c \u b2 \AN 

If site 1 has many such neighbors, then the transition rate to the continuum is multiplied by 
the number of neighbors. What happens is that the stae leaks first through an impurity and 
then into the continuum. Rare transitions from one low energy impurity state to another 
will therefore not lead to the formation of a well-defined band, and transport occurs through 
the extended states. 

The result is a collection of resonances together with a set of broadened continuum 
levels. As in the semiconductor case, the total number of continuum states is equal to the 
number of atomic orbitals in the sample; the number of resonances is equal to the number of 
impurities. The resulting total density of states is shown in Fig. [TJ]. We must conclude that 
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the low temperature behavior in this system is never dominated by interimpurity hopping. 
As a result, calculations using impurity averaging methods should lead to correct results. 
However, the noncrossing approximation may be dangerous in two dimensions, and the use 
of a symmetric density of states is not justified. 

I am very grateful to A. V. Balatsky and A. Chubukov for helpful discussions. This 
work was supported by the National Science Foundation through Grant DMR-92-14739 
and through the Materials Research Science and Engineering Center Program, Grant No. 
DMR-96-32527. 



VI. APPENDIX 



In this appendix I give details of the more lengthy integrations. Eq. |151 



is 



re ~\ N °r^fj f 1 ^ + ^3 + cos(0)(1-x 2 ) 1 / 2 t 1 ^ 

y(iu,rx) = — — dtp as e — — ,. j- — dx. (151) 

v ' An J J J-i Q 2 + s 2 + cos 2 0)(l - x 2 ) y ' 

The last term in the numerator gives zero on integration over 0. We again write the x 
integral as 

1 dxe- iQ + ST3 



i fi 2 + s 2 + cos 2 0)(l-x 2 ) 

ifl + sr 3 



dz e 



Vt 2 + S 2 + COS 2 (0) — Z 2 COS 2 (0) 



poo 

-i dy e 
Jo 



n 2 + s 2 + cos 2 (0) - (1 + iy) 2 cos 2 (0) 



+i H dv e iu (- l+iy) lVL + STs (152) 

Jo y fi 2 + s 2 + cos 2 (0)-(-l+i|/) 2 cos 2 (0)' K ' 

The poles lie on the real axis outside the contour, so the contour integral is zero. Thus we 

obtain a very simple result: 

+ sr 3 ) e iu - e~ iu 

1 = ~WT7 2 — ' (153) 

again because u » 1. Then 

Arckpr J J il 2 + s 2 

= ^0 [e ik * r (sgn (uj)t + r 3 ) - e~^(^n (lu)t - r 3 )] . (154) 
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This is Eq. [L34. 



Let us now carry out the integration for a polar gap in the general case. The Green's 
function is given by Eq. |137| : 



7V f Zn f f 1 

Q(r(smipx + cos ipz), iu) = / d(f> de dxe 

Air Jo J J-i 



N ° f 2n ^ f 1 ^iuJ" + er 3 + Aon(ax - by/1 - tf) 



+ e 2 + Ag(ax - by/1 - x 2 ) 



....... , .... , u t( t/ y .1. - ■ ,i 

(155) 



where a = cosip, b = cos sin?/', and u = (kp + e/vp)r. The coordinates for the integration 
have been rotated so that the polar axis defined by x = ±1 is along f. This is rewritten as: 

No f 27T f 

Q(r(smipx + costpz,iuj) = — / d(f> dsI(Q,a,b), (156) 

47rAo Jo J 

where 



7(0, a, b) = C ^u^ + STs + njax-bVT-^) 
K ! J-i Q 2 + s 2 + {ax - by/T^ 2 ) 2 K ' 

which may again be evaluated using the same contour as above. This procedure yields 
J iV + s z + [az — by/1 — z z \ z 



-i r dye*^ — + + Tl(a(1 + W) - / * = (1 + ^ 



O 2 + s 2 + [a(l + %) - 6^1 - (1 + i|/) 2 ; 

roo . . ^0 + sr 3 + r 1 (a(-l+^)-6 x /l-(-l + ^) 2 ) 

i I dye m (~ +iy ' - — - — - - (158) 

Jo Q 2 + s 2 + [o(-l + iy) - 6^1 - (-1 + iy) 2 } 2 

The poles lie at points determined by the equation 

ax - by/l -x 2 = ±i(tt 2 + s 2 ) l/2 . (159) 

There are four such roots: 

&(a 2 + b 2 + ,2 + ^2)1/2 m( ^2 + s2) l/2 

(a z + b z ) a z + b z 

Only the two with positive imaginary parts (x ++ and x |_) may lie in the contour, and then 

only if \Re x\ < 1. This is the case if \s\ < (a 4 /b 2 + a 2 — O 2 ) 1 / 2 = sq. The cuts produced 
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by the square roots can be chosen to be along (—00, —1) and (1, 00). They lie outside the 
contour. Performing the integrations leads to 

iVt + sr 3 + i\/Vl 2 + s 2 Ti 



I(n,a, b) = e(s - \s\)ne tux ++- 



. e ikFr in + sT 3 + aT l cisAor/vF 
k F r Q 2 + s 2 + a 2 



e-** r ig + sr 3 - an Aor/vF 
+% k F r n 2 + s 2 + a 2e • [ b } 

This must next be integrated over s. The first two integrals can be performed by the 
stationary phase approximation, taking u as a large parameter and expanding the argument 
of the exponential about its maximum. These expansions are valid if k F r sin if) » 1 and 
k F r cos if) » 1. Thus, the result is not valid near the equatior or the poles. Fortunately, 
we already have results in these regions. The second two integrals are standard contour 
integrations. This leads to the expression 



/ I(n,a,b) = 2m(sgn (u) + n)[ — — ^] 1/2 exp 

j— 00 au 



, r 2n \Q (a 2 + 6 2 ) nl/ o 

+2m(sgn (u) + n) — !— ^ ^ 1/2 exp 

aw 



i-ufr v 7 a 2 + o 2 + f2 2 — au\Q\ 
a 2 + b 2 

-iuby/ a 2 + b 2 + f2 2 — adQI 



a* 



+ 2 



2 



A; F r V^ 2 + a 2 V^ 2 + a' 

The final integration over the azimuthal angle is also simplified by the fact that u » 1, 
and the same method may be used. 
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N Q A f 2n f 
Q(r[smif}x + cos if)z,iu}) = / d(f> / dsI(Q,a,b) 

4:71 Jo J 



iN A /"^ 27r|fi|(cos 2 ^ + sin 2 ^cos 2 



(s^rt (fi) + n) / [— ^ — — 



exp[ 



2 Jo u cos ^ 

r iu sin if> cos — it|fi| cos^-. 



X 



cos 2 if; + sin 2 if; cos 2 - 



ziV(A W (n) + n ) / 2 ^N(cos 2 ^WyW0) ]1/2 x 



exp[ 



2 jo it cos ^ 

—ill sin if} cos — cos^ 
cos 2 ^ + sin 2 if} cos 2 - 



—iN Tre lkFr r c<j cos ^7~i , 

2/c F r X 2 + Agcos 2 ^ X 2 + A 2 cos 2 V 



e 



- ^ a; 2 +Ag cos 2 t/t /up 



iN 7ie lkFV UJ COS 

2A; F r X 2 + A 2 , cos 2 ^ V^ 2 + cos 2 V 



e 



-^/w 2 +A 2 cos 2 i/>r/i>F (163) 

The integrations then give 

g(r(sm?/>x + cosif}z,iu}) = ~~^T^( s 9 n + r i)[ A \ „„„„ ;J2 | „ ;J2 J x 

gikpr sinip—kpr\w\ | cos i/>|/Ao 

~2fc F r 

g— ifcj?r sin i/>— fei?r|u>| | cos i/i|/Aq 



2\u\ 


A 


cos?/' 2 


sin if> 2 




2|w| 




A 


cosif>\ 2 


sin if) 2 



to" M + n)[ A , 7' . /l2 ] 1/2 



—iN Tce lkFr , uj cos 

2k F r Juj 2 + A 2 , cos 2 if} Ju) 2 + A 2 cos 2 if} 



— y/ w 2 +Aq cos 2 ?/>r /t>ir 

iN ire' lkFr to cos t/>ti 



+ * ; = - r 3 x 



2 ^r ^2 + a 2 C os 2 if} ^u 2 + A 2 , cos 2 ip 

e - y/u) 2 + A 2 cos 2 Vr/ t> F _ (164) 



This leads immediately to Eq. |138 



The Green's function for motion away from the impurity along the x-axis in two dimen- 
sions is given by Eq. |140| : 



Q(iu rx) fde C dV ^ + ^3 + ^(2^-1) 



7T 
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Here v = cos(0), A = A (2v 2 — 1) and u = (/cp + ejvp)r. We may also write 

2iV n 



Q(icu,rx) = / dsI(Q,s), (166) 

7T J 



where 



I(Q s) = C dv ^ + ^3 + ^(2^-1) 

and Q = cj/A , s = e/A , and t = \/VL 2 + s 2 . As usual, we split the integral into three 
parts: 

s) = h(n, s) + i 2 (n, s ) + i 3 (n, s ), (168) 

with 

T (Q «\ - I dz ^ + ^3 + ri(2^ 2 -l) 

h^s)-f-j=== t2 + [2z2 _ iy2 e , (169) 
y 1 Jo . \ (\ _i_ ,-,,^2 t 2 + 2 1 + w 2 -l 2 ' v ; 



and 



l-(l + iy)2 t 2 + [2(l+iy)2-l]2 



W -) = i r , = ^ t/ 7 " 3 rt n J 2( ~\t ^ ^<-^>. (171) 



l0 ^l_(_l +iy )2 i2 + [2 (_l +iy )2_ 1] 

Turning first to ii, we locate the poles of the integrand at 



r + (2^ -1)^ = (172) 



or 



^ = i(l±it). (173) 

There are four poles, only two of which have positive imaginary parts and thus may lie in 
the contour. Let us call them z\ and z 2 . They satisfy 

N = N = M-^(i + * 2 ) 1/4 (174) 
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and z\ = — \z\e~ %e , z 2 = \z\e ld , with 9 = |tan _1 (t). If the real parts of these quantities are 
less than unity in absolute magnitude, then they lie in the contour. A bit of trigonometry 
shows that this occurs only if t < 2\[2. This condition will normally be fulfilled in our 
problem. The other two roots lie at z 3 = —\z\e' ld and z A = \z\e~ ld . 
The contour integral is now written as : 



dz 



iQ + sr 3 + n(2z 2 - 1) 



4 J y/1 - Z 2 (z - Zi)(z - Z 2 )(z - Z 3 )(z - Z A ) 

and the residue theorem gives 



(175) 



h(n, s) = ^2mG(V8-t) x 



Using 



in + st 3 + n(2z 2 - 1] 



+ 



iQ + sr 3 + r x {2z\ - 1) 



,JUZ2 



z\(z x - Z 2 )(Z 1 - Z 3 )(Z 1 - Z A ) z\(z 2 ~ Z 1 )(z 2 ~ Z 3 )(z 2 - Z A ) 



• (176) 



vT^ = _L(i + t2 )i/^ = 



z\e 



7 2 
z 2 



V2 



(l + t 2 f A e~ M = \z\e 



2z\-\ = -it, 



2zi 



1 = it, 



(177) 



this may be rewritten as 



h(n, s ) = ^e(V8 - t)\z\~ 4 x 



(in + st 3 - itr x )e iuz ^ 



e tf(_ e -tf - e i6 )(-e- ie + e w )(-e-* e - e- id ) 
(in + st 3 + itr^e™* 2 



e -i0( e ie _|_ e -i8)( e ie + e ie )(e w - e^ 6 ) 1 

|e(V8- t ) (T ^)x 

(in + st 3 - itr x )e iuzi (in + st 3 + itr x )e iu ^ . 
H-2cos0)(2i sin0) (-2) + (2cos#)(2)(2isin#) 



(178) 



Since \ sin 9 cos 9 sin 29 = sin tan 1 1 = t/\/l + t 2 , this becomes 
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h(Q,s) = ^0( v / 8-t)( 



x 



2 v /y tVT+¥' 

[(iVl + st 3 - itn)e iuZl + (iVl + st 3 + itn)e iuZ2 ]. (179) 



Integrating this: 



/7T roo I 
*/,(«, s) = - ^ - v TO)( 7 _ wm _ )x 

X 



(ifi + sr 3 — iV^ 2 + s 2 Ti) exp 
exp ^y/y/TTwTsZ-l) 

+ (itt + sr 3 + WW + s 2 n) exp ^V / v / TTfi 2 T^ 2 + l) : 



exp ( (^VVl + fi 2 + s 2 - 1 



(180) 

Since u » 1, we expand the argument of the exponential around its maximum at s = 0: 



exp(— vVl + ^ + s 2 - 1) « exp( — VVl + ^ 2 - 1) x 

exp[^(v / r + T^-l)- 1 / 2 (l + fi 2 )- 1 / 2 ], (181) 
8 

and evaluate the rest of the integrand at s = 0, so that t — > Q. This yields 

[(i s#n Q - iri) exp(— 

+(isqn Q + in) exp(^V Vl + ^ 2 + 1)]. (182) 
The other integrals are simpler: 



^y 2 - 2iy t 2 + [2(1 + iy) 2 - l] 2 
V ; t 2 + l Jo ^ 



1_i ikpr isAor/vF & + ST 3 + n 



and similarly 
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These expression are easily integrated over energy: 



r dsl 2 (n, s) = (1 - ^ t &-e-^^ + ^L+S* + Tl , (185) 

i-oo v ; 2 y k F r VTTTi? ' v ; 

and 

r dsHQ, s) = (1 + ^ .Ee-^e-^^ ^ ' *^ + ^ + n 
i-oo V ' ; 2 V^r Vl + fi 2 



47 



Combining Eqs. |165| , |182| , |185| , and |186| , we find 
2N 



Q{iu, rx) 



71 



\{i sgn n - m) ex P (— -\fVTTtt 2 + i) 



+ {isqn n + in) exp(y V Vl + ft 2 + 1)] 



, C 1 ~ c -ik F r c ^vr+^A r/v F in ~ iy/i + ^ 2 r 3 + n 
2 V^ r v/TTTF 



(187) 



From this Eq. |141| follows immediately. 



The Green's function for the direction along the node is more easily computed. We have, 
from Eq. [142], that 



r( - x + y 2N r f 1 dv iu + er 3 + A n*Vl - v 2 iuv 

Qiiu.r — =-) = / de / . — „ —. r^e , 188 

V ' y/2 ! 7T J J-i t/T^ u 2 + e 2 + 2A 2 v 2 (l-v 2 ) v ; 

where the part proportional to n vanishes by symmetry, and with the usual breakup: 

g( iuJ ,r^J) = -^ J ds(h + I 2 + h). (189) 



Now 



dz in + st 3 _^ iuz 



vT^I 2 n? + s 2 + 2z 2 (l -z 2 ) 



If l n + sr 3 + r 1 zVT^ (19Q) 



4 J (z - z x ){z - z 2 )(z - z 3 )(z - z 4 ) 



where the roots are 



Z2 = -Zi 

z 3 = -^{VTTT 2 + i)^ 

z A = ~z 3 . (191) 
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We shall also need: 



(192) 



Only z\ lies in the contour. Evaluating the integral: 



4 [Z! - z 2 ){z 1 - Z 3 )(Zi - Z A ) 



-iri in + sr 3 _ e%uzi 



2 z A {2 Zl ){zl - zj) 



Integrating over energy: 

iVt + sr 3 —u 



[ dsh = 71 T ds - | gT3 exp [^(Vl + ^ + g 2_ ^1/2] 

J Loo V^ 2 + sVi + n 2 + s 2 l Vz 1 J 

= n g eX p[^(v / rT^ - 1) 1/2 ] x 

/OO 7/S 2 
exp [^L- ( VT+Q? - 1 ) 1/2 VT+Q 2 } 

= 4tt 3 v^ 1/2 isgnQ ^/j-yn2 1 )i/4 c -fc f r( v W-D 1 / 2 /^ i 
fcpr (1 -|- ri^)^/"^ 

The other integrals are simpler, as usual. 

Jo ^1 - (1 + iy )2 Q 2 + s 2 + 2(1 + iy) 2 {l - (1 + iy)2) 



/•oo 



1-i e ikpr in + sr 3 



and 



2 s 2 + Q 2 



+ iy) 2 W + s 2 + 2{-l+ iy) 2 {l - (-1 + iy) 2 ) ' 



/•oo 



ie~ tu / -^e" ny 



_ u JVl + st 3 



2 s 2 + n 2 

The energy integrals are easily done: 



(194) 



o ^=2~w s 2 + n 2 

e ik F sA r/v F ^ ^g 5 ^ 



o ^Ha) s 2 + Q 2 

1+i in±STs p - ikFS A r/v F _ (1Q6) 
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/ 



dsl 2 = (1 ^ ^= (i sgn + *T 3 )e- Mr/ ^, 
2 V^r 



(197) 



and 



/ dsh = ^fi 1 & -^= (< s»n w - zr 3 )e-l-l^. 



2 v^f^ 

From Eqs. |189j , |194| , |197| , and |198| , the Green's function is 



(198) 



x + y, 
V2 ' 



2Na 



X 



7T 



+ (l_-z)vr ( . ^ ^ + ^-Mr/,, 



2 y/kpr 
(1 + z)tt e-^ r 



2 v^fT 
From this, Eq. |143| follows immediately. 



(z s^cj-ir 3 )e- |Lj|r/,;F ]. 
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The general direction is given by the equation: 
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The two integrals come from the regions = (0, tt) and = (71*, 2ir). Then we find 

7T iQ + ST3 + i|i|rx 
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2 |t|vTT^ 

where the result is obtained by the same method as that in Eq. |193| . The energy integral is: 
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Similarly, from Eqs. |1971 and |198| , we have 
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Finally, let us define 
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G-(iuo, ip) = G(iu, r cos(— — ip)x + sin(— — ip)y). 



Collecting results, we have 
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This is Eq. |45. 
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FIGURES 

FIG. 1. Density of states for a model semiconductor with a band structure which is symmetric 
about the gap midpoint. This system contains a single impurity which is a unitary scatterer. This 
results in a midgap state. 

FIG. 2. Density of states for a model semiconductor with a band structure which is not sym- 
metric about the gap midpoint. This system contains a single impurity which is a unitary scatterer. 
This results in a bound state which is not at midgap. This is the generic case. 

FIG. 3. Density of states for an artifical semiconductor with nearest-neighbor tight-binding 
dispersion and a gap. This density of states has an asymmetry factor of 0.25 t / unit cell. The 
asymmetry factor is defined in the text. 

FIG. 4. Density of states for a semiconductor with many impurities using the usual impurity 
averaging procedure. The impurity density is above the critical value, so the gap has filled in. The 
even density of states is for a semiconductor with a band structure which is symmetric about the 
gap midpoint. The parameters, referring to Figure 1, are ei/A = 5 and ni mp V 2 No/A = 0.8. The 
asymmetric denstiy of states is for an asymmetric density of states. The parameters, referring to 
Fig. 2, are e £ /A = 5, e u /A = 8, and n imp V 2 N /A = 0.8. 

FIG. 5. The density of states for a model semiconductor with many impurities. The calculation 
combines the results of the T-matrix and the impurity averaging. The concentration is subcritical, 
so the gap is not completely closed. The bound states are shown in black. Note that the total 
number of bound states is proportional to the number of impurities, whereas the number of states in 
a fixed energy range around the chemical potential is proportional to the total number of orbitals. 
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FIG. 6. The density of states for a model semiconductor with many impurities. The calculation 
combines the results of the T-matrix and the impurity averaging. The concentration is supercritical, 
so the gap is closed. The bound states, shown in black, become resonances. Note that the total 
number of states in the resonance peak is proportional to the number of impurities, whereas the 
number of states in a fixed energy range around the chemical potential is proportional to the total 
number of orbitals. 

FIG. 7. Contour for the angular integrals involved in calculating the impurity wavefunction in 
real space. 

FIG. 8. Diagram (a), which contains no crossed impurity lines, is counted in the usual calcula- 
tions of transport properties in unconventional superconductors. Diagram (b), with crossed lines, 
is usually neglected. 

FIG. 9. The diagrams of the previous figure in momentum space in the normal state. The 
justification for the neglect of (b) is that one of the momenta is off the Fermi surface unless an 
additional constraint is applied. 

FIG. 10. The same diagrams in momentum space, but now in the superconducting state at low 
temperatures. All momenta must be near the nodes, which are situated on the diagonals. The 
justification for the neglect of (b) no longer applies. 

FIG. 11. The density of states for a d-wave superconductor with many impurities. The calcu- 
lation combines the results of the T-matrix and the impurity averaging. The broadening of the 
resonance is shown. Note that the total number of states in the resonance peak is proportional to 
the number of impurities, whereas the number of states in a fixed energy range around the chemical 
potential is proportional to the total number of orbitals. 
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